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Abstract 

A method is presented for numerical implementation of the ex- 
tended TDHF theory in which two-body correlations beyond the mean- 
field approximation are incorporated in the form of a quantal collision 
term. The method is tested in a model problem in which the exact 
solution can be obtained numerically. Whereas the usual TDHF fails 
to reproduce the long time evolution, a very good agreement is found 
between the extended TDHF and the exact solution. 
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1 Introduction 



The study of nuclei in out-off equilibrium configurations in many-body quan- 
tum theory is a very complicated problem. However, due to their long mean- 
free path inside the nucleus and to their large amplitude zero point motion 
associated to a strong derealization of their wave-packet, nucleons could be 
often regarded as independent fermions moving within a mean potential. In- 
deed, such theory, called Time-Dependent Hartree-Fock theory (TDHF), first 
applied to nuclear dynamics some twenty years ago by Bonche et al 0, is 
able to reproduce quantitatively many properties of Heavy- Ion collisions p|. 
However, in this theory, nucleons interact only through the mean-field and 
collisions between particles are neglected. It is difficult to believe that corre- 
lations induced by collisions do not play also an important role when highly 
dissipative processes are involved. For example, thermalization of single- 
particle degrees of freedom || and damping of collective motion 0, ||] could 
only be understood by the introduction of two-body correction. Furthermore, 
usual mean-field theory is not able to reproduce the large experimental width 
of mass produced during Heavy-Ion reactions at Intermediate energies. 

In the extended Time-Dependent Hartree-Fock (ETDHF) theory, the de- 
scription is improved beyond the mean-field approximation by incorporating 
correlations in the form of a quantal collision term||. Due to the numeri- 
cal complexity, the applications of this theory on realistic situations remain 
a difficult problem, and therefore only a few approximate calculations have 
been performed so far0. In this paper, we describe the basic features of the 
ETDHF approach, propose a method for obtaining numerical solutions of 
the theory, and present an application of the method to an exactly solvable 
model. 



2 Extended Time-Dependent Hartree-Fock The- 
ory 

The exact description of a quantum system is contained in the many-body 
density operator D(t) and its dynamical evolution is given by the Liouville- 
von Neumann equation 



.,dD 

in—— 

dt 



H, D 



(1) 
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Since the system has a a priori a large number of degrees of freedom, the 
practical resolution of equation (JJ) is often impossible. The basic assumption 
of mean-field theory is to suppose that only one-body observables are impor- 
tant in the description of the system (neglecting , by the way, all correlations 
of order greater or equal to two) . Usual mean-field theory neglect completely 
two-body dynamics and values of two-body observables are obtained assum- 
ing that they are equal to the least biased value when one-body observables 
are known 0. In this approaches, all relevant information is contained in the 
one-body density operator p and the dynamical evolution is replaced by|§ 

*f=M (2) 

where h is the mean field associated to H. This approximate theory is al- 
ready a many-body theory for out-off equilibrium finite system. Note that, 
even if two-body dynamics are not explicitly included, non trivial correla- 
tions between single particle states exists due to the reorganization of the 
mean-field during the time evolution. 

However, this two-body correlation is not sufficient to reproduce the many 
facet of nuclei desexcitation. In particular, no dissipation, characteristic 
of the irreversible flow from collective motion to single particle degrees of 
freedom || and expected to be responsible of the thermalization of nuclei, 
is included. The next step towards extended mean-field is to include in- 
formation about two-body correlation dynamics. This could be obtained in 
particular by a truncation of the BBGKY hierachy||10|| to the two first equa- 



tions. This approaches assumed that TDHF is already a good approximation 
of the dynamics and two-body correlations are supposed only to act as a small 
perturbation added on top of the mean-field. This leads to a correction factor 
in the evolution of the one-body density operator[|[|ll[| 



dt 



h,p +K(p) (3) 



with 



K{p) = -ij dt'tr 2 [V 12 , U 12 (t,t>) PlP2 V 12 (1 - p x ) (1 - p 2 ) U l2 (t',t) 

-U l2 (t, t') (1 - Pl ) (1 - p 2 ) Vn PlP2 U 12 (f , £)] (4) 



In this expression, usually, an additional term exist due to the propagation of initial 
correlations. In this paper, we suppose that, initially, we have uncorrelated states. 
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Taken same notations as in ref. [TT], the label 1 and 2 refers to first and second 
particles, for example 



(ij\pip2\kl) = (i \p\ k) (j \p\ I) (5) 

Uu is the propagator of two independent particles associated to the mean- 
field: 

U 12 = Ux®U 2 (6) 

with 

Ui(t, = T (exp (- % - jf h(p ( a ))<fc)) (7) 

and tr?, is the partial trace taken over the second particle whereas V\2 includes 
the antisymmetrization. 

Eq. (f|) is the starting point of most of the theories that goes beyond 
TDHF. For example, all semi-classical applications which include collision 
effects [0] (namely BUU approaches) can be derived from it, by developing 
the Wigner transform of (^) up to first order in h. Over the past decade, 
significant progresses have been made on the description of nuclear dynamics 
within semi-classical transport theories and extension of mean-field seems to 
be a promising tool for describing dynamics. However, nuclei are quantum 
objects even under extreme conditions and important physical component 
are missing in a semi-classical treatment. Some attempts in order to treat 
the collision term in a quantum picture already exists]?], [L4]]. Nevertheless, 
due to the complicated expression of K(p), theses investigations have been 
carried out under strong approximations. 

Here, we describe a procedure for treating the collision term ([|) with 
minimum bias and which may be adapted for solving the ETDHF equation 
in realistic situations under reasonable numerical approximations. 

3 Re-examination and numerical algorithm 

In a single particle bases, the one-body density operator reads 

P = £iA)n Ai v<A'| (8) 

A, A' 
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and the evolution equation take the form|TJ 



ih4 i \x) = h{ P )\x) 



(9) 



with 



iV = E W\V l2 \af3) f dt' (a^lVi^; 

a,a',P,/3' ,7,5,5' 00 

(n 7 A' n «'(^«' _ *w)(<W - n W /) (10) 

-1o«' ™/3/3' (5 7 A' - n lX > ) (<W - «55' ) ) (11) 

This expression is the most general one could extract from (H). We know 



from work of [14], that approximations in (|TTD leads to important change in 
physical results. In particular, the time integration is often replaced by a 
conservation of energy based on a Markovian approximation of the collision 
but we know that this memory effect is important since it creates the non- 
trivial coupling between collective and single particle motions. Furthermore, 
we are obligated to consider the non-diagonal part of ( |TTD which will give the 
coherent evolution of single-particle states. 

A direct resolution of (|TT]) needs a large numerical effort. However, in the 
following, we will see that it could be considerably simplified without loosing 
its generality. 

3.1 Coarse-Graining in time 

Supposing that we know the one-body density at time t: 

p = Y.\®i) n i(®i\ ( 12 ) 

i 

From information theory, it means that we have access to the least biased 
description of our system when only one-body observables are known. We 
can have a geometrical picture of the information reduction. One-body space 
could be represented by a manifold in a bigger space (see fig. [TJ). Suppose 
that we start from one point on this manifold. The exact dynamical evolution 
could be represented schematically by a trajectory in this space. Due to the 
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presence of two-body correlations, the trajectory will not remain in the one- 
body space. Nevertheless, the exact trajectory is associated to a path in 
the one-body space which reproduce the one-body dynamics of the exact 
evolution. The goal of using mean-field theory is to predict this trajectory. 
In the one-body space, the TDHF theory corresponds also to a trajectory. 
However, due to the neglected two-body correlations dynamics, usual mean- 
field is not able to reproduce the long-time evolution of the system. This is 
symbolized in fig. [1], by a strong departure of the TDHF theory from the 
best one-body trajectory. The goal of ETDHF is to correct the mean-field 
evolution in order to be a much better approximation of the best trajectory. 




Figure 1: Geometric picture of TDHF and EXACT evolution. The Initial 
condition is supposed to be contained in the one-body space. The exact 
trajectory is symbolised by a path in the total space. In particular, due to 
the presence of two-body correlations, this trajectory do not remain in the 
one-body manifold. The best associated trajectory in the one-body space is 
represented in dashed-line whereas the TDHF evolution is represented in solid 
line. The neglected two-body correlation dynamics in TDHF, implies that 
this theory do not corresponds to the best trajectory for long time evolution. 

However, some remarks are in order: 
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• Whereas TDHF fails to reproduce the long time evolution due to accu- 
mulated errors on two-body observables, one-body dynamics remains 
dominant. It implies that it exists a typical macroscopic time interval 
At during which mean-field is a good approximation of the dynamics. 
This time interval is of the order of the time between two collisions. 
Note that, a basic assumption of ETDHF and TDHF is that this time 
interval is large compared to the typical time-scale of mean-field evo- 
lution. 

• For time greater or equal to At, two-body correlations appear as a small 
correction to the mean-field and could be treated in perturbation. 

It enable us to coarse-grain the time evolution for the extended mean- 
field dynamics. Instead of solving equation (^) in time, we will divide the 
evolution into two step: 

• The Hartree-Fock states will be evolved through the usual TDHF equa- 
tion between t and t + At 

■ zh§- t \<i> t ) = h(p)\^) 

(13) 

{ im = a 

• We then use the perturbation theory in order to correct the one-body 
density by the error accumulated during At associated with (|9]). We 
then find the new density operator //(see fig.§) 

+ Ap 

ft 

Integrated (14) 

Effect 
of collision 

The new one-body density is then diagonalize and the procedure is 
iterated (see fig. §). 



p'(t + At) = p(t + At) 
ft 

th§ = [h,p] 
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Figure 2: Schematic representation of ETDHF evolution within the coarse- 
graining time framework. The two steps procedure is represented. In solid 
line, the usual TDHF is first performed. After a time At, a correction is 
applied (small dashed line) in order to account for two-body correlations 
propagation. 

3.2 Determination of a Basis 

In order to express different operator matrices elements, we need a priori 
a complete basis of the single-particle space. However, this could rarely be 
obtained in a problem with a huge number of degrees of freedom. On the 
other hand, if we suppose to have such a basis, only few of the states will 
be necessary in order to have the main information about the system (the 
rest corresponding to configuration which are not accessible). In a previous 
application 0], the author proposes to include only a part of the Hartree- 
Fock bases (including all occupied states ("hole states") and few unoccupied 
states ("particle states")) and to follow them with time. However, the system 
could reach many different configurations and single-particle states that were 
not included at the initial time could become important. It is thus difficult 
to believe that one can choose the basis at the initial time. Conversely, at 
each coarse-grained time step At, we will construct a bases that contains the 
important information about different configurations accessible to the system 
between a given time t and t + At. According to this procedure, 

• At a given time t, we suppose to have the Hartree-Fock basis |$j(t)). 
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Exact evolution 




Correction 



Figure 3: With the same convention as in fig.l, the schematic evolution 
through the numerical resolution of ETDHF is represented. Each At, a 
correction is applied to the mean-field evolution. The ETDHF procedure is 
expected to follow the best trajectory in the one-body manifold. 
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We only perform the evolution of occupied states between t and t + At. 

• At time t + At, we complete the hole states ( \&i(t + At))) by a limited 
ensemble of particle states( \pj(t + At))). This states are constructed 
in order to be a good approximation of eigenstates of h(t + At) with 
an energy smaller than a given energy Euax- in the following, we will 
discussed in detail the construction of the basis and the choice of Smox- 

=3- Between t and t + At, two single-particle states (1-2) could collide. 
As a consequences, new states could be populated. In a simplified 
vision of collision^], we can assume that the transition probability 
from (1-2) to two new states (3-4) is important when 

E\ + 82 ~ S3 + 64 

The maximum accessible energy Euax is obtained when 1 and 2 
are the hole state with highest energy Euax and 3 or 4 is the hole 
state with lowest energy Euin- This give the natural truncation 
energy for the basis: 

£Max — 2Emclx — Emiti 

=>- In order to find states with low energy, we apply an imaginary time 
method to Hartree-Fock states. This method, explained in AP- 
PENDIX A, provides finally a truncated basis: {\§i{t + At)) <g> \pj(t + At))} 
where all different matrices could be expressed. This basis will be 
called "instantaneous basis" and noted generically {|A)} in the 
following. 

• In order to include the memory effect in the collision terms. The in- 
stantaneous basis is evoluted backward self-consistently. 



4 Re-organization of occupation numbers and 
single-particle wave-functions. 

The collision term introduce a small correction to the density at each macro- 
scopic time step At, see expression (14). In the following, we will treat Ap in 

2 Here, the simplification is contained in the fact that, in general, states (1-2) are not 
eigenstates of h which is the condition of having an exact conservation of energy. 
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first order perturbation theory. In this case, the new density operator could 
be expressed as: 

P' = EI$aK<^I (15) 

i 

where new states are obtained through the following procedure 

• New occupation numbers: In first order perturbation, only diagonal 
elements of Ap are necessary to calculate the new eigenvalues n' x of p'. 
The time evolution of this diagonal elements between t and t + At is 
given by a master equation (see APPENDIX B) 



dt 

with initial condition 



f/ " A (l-n x )Wt-n x W; (16) 



n' x {t) = n x (t) (17) 



During t and t + At the gain and loss term, respectively W x and 
W x , are considered constant. In this case, the equation ([UJ) is exactly 
solvable and reads: 



n x >(t + At) = n x (t) exp (-At (w£ + W 



A 

+ wfw t 1 - exp (~ At W + (18) 

A A 

Reorganization of states: The new states are given in perturbation 
theory by 

m + At)) = |A> + £ |V) (V| Ap |A> (19) 

It is important to note, that, not only previously occupied states but 
also part of unoccupied states are included at each macroscopic time- 



step in the dynamics. In a previous application [14], Tohayama argued 
that non-diagonal elements in the density matrices are important. In 
our approach, they are explicitly included since the first order pertur- 
bation is equivalent to a approximate diagonalization of p' . However, 
here we do not fixe the basis at the initial time and the system will 
dynamically choose which configuration will be accessible. 
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In this section, we have described the scheme we will use in order to 
extend TDHF. In the following, we will apply it to a model that could be 
exactly solved numerically. 



5 Application 

We have considered two distinguishable particles in a one-dimensional an- 
harmonic potential. The total Hamiltonian reads 

^ = E(^ + ^)+E^ (2°) 

where the one-body external field is 



Ui = \kx\ + \k'xt (21) 



and the non-local two-body interaction Vij is taken as 

^ = 7 SUexp-<fi--^ (22) 

This problem, which seems simple in appearance contains already different 
features of many-body physics and is a good benchmark for ETDHF since it 
could be compared with the exact solution. 



5.1 Initial conditions and evolution 

We consider the system initially heated and constrained. The initial two- 
body density operator is written as a statistical equilibrium 

£,i " i = i exp -^( J? - A °) (23) 

where Q is a one-body operator and XQ is the initial constraining field. At 
initial time, the constraint is relaxed and the system evolve. In fig. |], we 
show the one-body external field part of H in r-space (solid line) and the 
equivalent constrained field H — AQ(dashed line). After relaxation of the 
constraint, three different evolutions are compared: 
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Figure 4: One-body part of the external field in r-space, with the initial 
constraint (dashed line) and without the initial constraint (solid line). 
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EXACT evolution: For the exact evolution, we solve exactly the von 
Neumann equation for the two-body density operator: 



.+ dD 
in—— 
dt 



H,D 



(24) 



TDHF: The equation of motion is 



»*^ = [M (25) 
Note that occupation numbers remains unchanged during the evolution 

dm , . 

— - = 26 
dt v ; 

• ETDHF: For the ETDHF, we apply the procedure described above. 
In this case, occupation numbers evolve in time and new states are 
mixed to initial one during the dynamics. 

For both TDHF and ETDHF, the initial density is the one-body that 
corresponds to the exact one 

p ini = tr 2 D ini (27) 



5.2 Results and Discussion 

5.2.1 Prediction of one-body dynamics: 

In fig. [|, we have represented the dynamical evolution of the diagonal part 
of the one-body density in r-space. At initial time, the one-body density has 
been imposed to be the same for all calculations. We see that usual TDHF 
(dashed line) calculation is a good approximation of the exact one-body 
dynamics (circles) for time smaller than 100 fm/c. However, for longer time, 
we observe strong differences between TDHF prediction and the expected 
result. On the same graphic, the equivalent evolution is reported for ETDHF 
simulations (solid curve). The ETDHF prediction of the one-body density 
operator follow closely the exact one even for long time evolution. 

Since the one-body density matrices govern the evolution of all one-body 
observables, it means that the extension of mean-field enable us to predict 
one-body dynamics with a good accuracy. Indeed, if we look to one particular 
one-body observable, for example the center of mass motion (denoted (X), 
fig. we observe a very good agreement between exact calculation and 
ETDHF prediction whereas TDHF fails to reproduce the long-time dynamics. 
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Figure 5: Evolution of the one-body density operator in r-space. Three 
different simulations are displayed: Exact (small circle), TDHF (dashed line) 
and ETDHF (solid line). 
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200 400 600 



t(fm/c) 



Figure 6: Center of mass motion (< X >)X of the two-particle in interaction. 
The three different simulations are represented: Exact (small circle), TDHF 
(dashed line) and ETDHF (solid line). 
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5.2.2 Variation of occupation numbers 

In fig. [7], we have plotted the variation of occupation numbers (i.e. the eigen- 
values of p) in the exact (circles) and ETDHF treatment (dashed line). Due 
to the presence of a residual two-body interaction, we observe a reorganiza- 
tion of occupation numbers in the exact evolution. The numerical procedure 
we have developed in order to extend the mean-field dynamics is able to 
reproduce this complicated behavior. This demonstrate in particular that 
two-body correlation dynamics is properly taken into account through the 
method we have used. It is important to note that reorganization of occupa- 
tion numbers is accompanied by an evolution of the relevant single-particle 
states and that no bias exists in our approach, in the choice of this relevant 
states. This is one of the major improvement contained in our procedure. 

5.2.3 Conclusion and outlook 

In this article, we have discussed the possibility of extending mean-field in 
quantum dynamics in order to include two-body correlations. Such theory, is 
numerically much more complicated to apply than usual TDHF simulations. 
As a result, only few applications have been carried out up to now. These ap- 
plications often rely on important simplifications made in the collision term. 
We propose here, a procedure to solve numerically this problem. In our ap- 
proach, no additional bias is introduced by the simulation and the collision 
effect is treated in all its generality. The method is then tested on a simple 
model that could be exactly integrated numerically. This model of two parti- 
cle interacting one with the other already contains many facet of many-body 
problems. In particular, the two-body interaction leads to a strong reorga- 
nization of occupation numbers during the dynamics and affect considerably 
the long-time evolution. As expected, usual mean-field dynamics is not able 
to account accurately for this complicated behavior. On opposite, the in- 
clusion of two-body correlation dynamics into the mean evolution improve 
considerably the prediction of one-body observables. 

The ETDHF seems to be a very promising tool for the description of 
nuclei under extreme conditions. In forthcoming work, we will apply the 
numerical procedure to realistic nuclei. The inclusion of two-body correla- 
tions will considerably improve our knowledge about the nuclear dynamics. 
In particular, this many-body theory includes the possibility for collective 
motion to dissipate energy through an irreversible flow to single-particle de- 
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Figure 7: Occupation numbers evolution. Circle represents the exact evolu- 
tion of occupation numbers and dashed line represents evolution predicted by 
ETDHF. The TDHF simulation (not represented) predicts fixed occupation 
numbers. 
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grees of freedom. This will certainly ameliorate the prediction of damping of 
giant resonances. On the other hand, this extended theory is needed for the 
understanding of thermalization of nuclei. 

6 APPENDIX A: Determination of the in- 
stantaneous basis. 

The instantaneous basis is constructed in order to be a good approximation 
of part of the eigenvectors of the mean-field Hamiltonian h at time t + At. 
Noting |£ Q ) this eigenvectors, we have 

h\p(t + At)]\S a )=e a \S a ) (28) 

Determining all \£ a ) is not possible in general since it imply the inversion 
of huge matrices. However, the |£ Q )'s form a complete basis of the single 
particle space. In particular, Hartree-Fock states could be expressed in this 
basis 

l*i> = E<2l&> (29) 

a 

Using the imaginary time method, consists in applying the operator exp (—f3h) 
to the state: 

|$;> = exp (-/?/*) |$,> 
where (3 is a real number. The interest of such operator leads in the following: 

(exp (-(3h)) n \<5> t ) = exp (-(3ne ) £ C l a exp {-f3n (e a - e )) |Q (30) 

—a 

which shows that 

(exp (-/?/!))" |$ 4 >?=^|eo> (31) 

We thus see that, application of the imaginary-time operator to any state 
converge towards the lowest state in energy. More generally, this operator 
remove the space pertained by states with high energy on profit of low energy 
states. Using this property, new states are constructed with following steps: 
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Application of exp (—/3h): from each Hartree-Fock state, a series of 
state is constructed 3 : 



^ {0) > =1*0 

(32) 

^)) = 7V (n) {exp (-(3h) - exp (-/%) 

The number of states per Hartree-Fock states will determine the pre- 
cision of the evaluation of the collision term. In practice, we observe 
that only few extra states are necessary to be included between two 
times separated by At. 

All states are grouped and orthonormalized. Note that the Hartree- 
Fock basis remains unchanged. 

Instead of considering directly this basis, we diagonalize the projected 
Hamiltonian: 

(1-P)h[p}(l-P)\p j ) = s\p j ) (33) 

with 

p = Y^mm (34) 

i 

Finally states with energy e greater than the energy Euax are rejected. 
Note that, all created new states belongs to the kernel of pit + At) 
which will simplify considerably the expression of the collision term. 



7 APPENDIX B: Expression of transition el- 
ements. 

In our instantaneous basis, the density is diagonal, the coefficients of eq. |TT| 
reads 

iV = E (M\Vu\aP) A \ t 

3 Note that, the use of a field hi = h + Sh where 5h is a small external stochastic field 
improve this method. 
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t 

dt' (n x ,n s (l - n a )(l - rip) - n a n p {l - n y )(l - n 5 )) 



(af3\V l2 \X'S) 



A It' 



In order to express the integral in time, we applied a backward mean-field 
evolution (where we have neglect the correlation part). In this case, oc- 
cupation numbers are constant, however, a finite life-time r (equivalent to 
the Landau beating time) has been added in order to properly calculate the 
integral (we have fixed r in order to not biased the integral value). 

• Diagonal elements: The equation of motion is written as equation 
(16) where 

4 

W~ = -2 n oJifi{l - n s )Real{ (\5\V 12 \af3) A \ t 



J' dt' exp ( t —t ) (a/3 1 V12 1 XS) 



A \t' 



W+ = pEi 1 - " n p )n s Real { (X8\V 12 \af3) A \ t 

a,0,6 



/^rffexp ^— — j (a(3\V 12 \X5) A \ t j 

We have checked that the lost term and the gain term are almost con- 
stant during each time interval At 

• Non-diagonal elements: The non-diagonal elements are necessary in 
order to express the mixing of states. The integrated effect of collision 
is simply taken as 

(A|Ap|A') = ^(iV+i^ A ) (35) 
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